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ABSTRACT 

The passivity theorem may be used to design robust controllers for structures 
with positive transfer functions. This paper extends this result to more general configu- 
rations using dissipative system theory. A stability theorem for robust, model-independent 
controllers of structures which lack collocated rate sensors and actuators is given. The theory 
is illustrated for non-square systems and systems with displacement sensors. 
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1 Introduction 


In this paper we are concerned with control of systems with second-order dynamics modelled 
by the following system of ordinary differential equations. 

Mx(t ) + Dx(t ) + Kx(t ) = Fu(t) (1) 

y(t) = C d x{t) + C v x(t) + C a x(t) (lb) 

For structures, this is an approximation to an (infinite-dimensional) partial differential equa- 
tion model and this system of equations may be of very high order, although the structural 
matrices M , D and K are positive definite and typically sparse and symmetric. 

Control theory for time-invariant linear systems which are described by first-order 
dynamic equations has been well established for decades, and many control software tools 
are available today for systems written in first-order forms. For applications, engineers can 
simply convert whatever models they have to the first-order forms and then use the existing 
tools to design the controllers. 

For second-order systems, transforming to first-order form not only increases the 
dimension of the problem, but also destroys the sparsity of the structural matrices. Not only 
physical insight, but computational efficiency is often lost in conversion to first-order form. 
Existing control analysis and design software may not be able to handle such a large system. 
For example, solving a 1000-by-1000 Riccati equation is practically impossible with today’s 
numerical techniques. 

There are basically two ways to address the controller design problem for a large 
scale system. One way is to minimize the dimension of the system model, through some 
model reduction technique. The reduced model is used in the controller design and some 
robust design methodology is used so that the controller stabilizes the original high order 
model. This is the approach behind most H ^ design for structures [6]. 

Another way is to design a controller which is independent of uncertainties in 
the system model. The advantage to this approach is that, unlike the first approach, the 
accuracy of the original high-order model (1) is not assumed, nor is an accurate system 
identification required. This is very appealing, since at the current time, no accurate model 
of structural damping exists, and the existing models of stiffness are only accurate at low 
frequencies. Furthermore, damping is difficult to identify experimentally. 

Structures with collocated rate sensors and force actuators (C' v = F ) are passive, 
i.e. have positive transfer functions [2]. The passivity theorem [1, 3] implies that any positive 
controller will lead to a stable closed loop system. This fact can be used to design control 
systems which remain stable despite large modelling errors. Recently, several researchers 
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have incorporated this approach with other robustness theorems in robust controller design 
for systems which have collocated rate sensors and force actuators but whose transfer func- 
tions may be non-positive due to computation delays, actuator dynamics etc. [7]. However, 
the restriction of collocated sensors and actuators is a stringent one, and excludes many 
applications. In particular, systems with displacement sensors and systems with a different 
number of inputs and outputs cannot be handled with this approach. 

Since altering the configuration of sensors and actuators does not alter the fact 
that structures with positive damping dissipate energy, the question arises as to whether 
the passivity theorem can be generalized. For instance, when a mass-spring-dashpot is 
attached to any mechanical system, including flexible space structures, the damping of the 
system is almost always augmented regardless of the system size. The parameters of the 
mass-spring-dashpot are arbitrary, model-independent and thus insensitive to the system 
uncertainties. Knowledge of the system model may be used in adjusting the parameters in 
order to satisfy system performance requirements. However, changes or errors in the system 
will not destabilize the system because it is an energy- dissipative device. The question 
arises as to if there are any feedback controller designs using sensors and actuators which 
behave like the mass-spring-dashpot device. The answer is positive as shown in [4], wherein a 
conventional engineering approach was used to design active controllers which behave like the 
mass-spring-dashpot device. Although the conventional approach provides physical insights, 
it lacks the theoretical basis necessary for application to more complicated systems such as 
nonlinear flexible robots. This paper is motivated by the desire to provide a rigorous theory 
to support this argument based on physical intuition. 

In the next section we summarize results on dissipative system theory obtained 
mainly by Willems [8, 9]. A stability theorem which generalizes the passivity and small 
gain theorems is given. In order to facilitate the application of these ideas to non-linear 
models of structural dynamics, no assumption of linearity is made. In subsequent sections 
we apply this theory to derive conditions for robust, model insensitive control of (1) non- 
square systems with rate sensors (2) systems with displacement sensors. As an example, we 
present a rigorous proof of the stability of a mass-spring-dashpot type controller. 

2 Dissipative System Theory 

Intuitively, a dissipative system is one which dissipates energy. In other words, the energy 
stored in such a system will be less than the energy supplied to it. This concept can be 
written formally in terms of a storage Junction which generalizes the concept of energy. The 
theory in this section is taken from [5, 8, 9]. 
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Definition 2.1: A dynamical system E is defined as follows: 


0 ) 


IA is the input space and consists of a class of 17-valued functions on the positive real 
line. 


(ii) y is the output space and consists of a class of Y- valued functions on the positive real 
line. 


(iii) Both Y and U are finite-dimensional inner product spaces and both y and U are closed 
under the shift operator i.e. if u[t ) € U then u{t + T) € IA. 

(iv) Define := {(^j^i) € RxR\t 2 > t\}. The state transition function (f > : R 2 xXxU i— » 
X and defines the state through the relation x(t 2 ) — </>{t 2 ,ti, x (ti), u ) This function 
satisfies the usual axioms for autonomous dynamical systems: 


(a) Consistency of initial condition: 4>{t, t, xo, u) = x 0 

(b) Semigroup property: <f>(t 2 , t x , t 0 , x 0 , u), u) — <f{t 2 , t 0 , x 0 ,u) 

(c) Causality: ui(t) = u 2 (t) for t 0 < t < U implies = <f>{t,to,x 0 ,u 2 ) for 

< i < *1 . 

(d) Time Invariance: <j>(ti + T,i 2 + T,x 0 ,uf) = <f>(t 1} t 2l x 0 ,u 2 ) for all T > 0,t 2 > 
h,u 2 (t) = ui(t + T ). 

(v) The function r : X x U m > Y is the output function: j/(t) = r(x(t), u(t)). 

In view of the system time-invariance (4d) we will henceforth use to = 0. 

Definition 2.2: A dynamical system is said to be dissipative if there exists a 
nonnegative function S:Xh R + , called the storage function and a supply rate w : U X Y h- » 
R such that , for all fi >0, 

S(x 0 )+[ w(u,y)dt > S(x 1 ) (2) 

Jo 

where u EU,x i = <f>(ti,0,x O) u) and y = r(x,u). 

Note that the storage function is, in general, not unique. The following theorem 
shows that if an appropriate supply function is found, it is not necessary to actually define 
a storage function. 


Theorem 2.1 [8] If there exists a supply function w(u,y ) for E such that 

rh 

/ w(u,y)dt > 0 
Jo 

for all t\ > 0 then there exists a storage function S such that (2) is satisfied and so E is a 
dissipative system. 
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Moylan and Hill [5] define a useful class of dissipative systems . 

Definition 2.3: [5] Given P, Q, R of appropriate dimension with Q, R symmetric, 
define the supply rate 

w(u,y) = (y,Py) + (y,Qu) + (u,Ru ) (3) 

A system is (P, Q, P)-dissipative if 


r* i 

/ w(u,y) > 0 Vti > 0. 

JO 


Passivity, or positivity, can now be seen to be the special case where P = P = 0 
and Q = I. Strict passivity is obtained by defining R = —S^I and Q = I. Similiarly, we 
obtain input-output stable systems by choosing P = —I and R = k 2 I where k is the gain. 

Suppose G is a given system, for which we wish to design a controller H, arranged 
as shown in Figure 1. 

The feedback system, or alternatively the pair ( G , H ), is said to be externally 
stable if ux,u 2 €E £ 2 ( 0 , 00 , U ) implys yi,y 2 6 P 2 (0> oo,!^)* and there is a maximum ratio, the 
L 2 gain, between the norm of the input and the norm of the output. 

The following theorem provides a simple test for external stability of intercon- 
nected (P, (5,P)-dissipative systems . 

Theorem 2.2 [5][Theorem 1] Consider systems G and H connected in the familar feedback 
configuration shown in Figure 1, and assume that H is (P 2 , Q 2 , R 2 )-dissipative and that G 
is (P 1 ,Q 1 ,Ri) -dissipative. The closed loop system is externally stable if 

a ._ + P 2 Qi - Qi / 4 \ 

Q' 2 -Qi R2 + P1 w 


is negative-definite. □ 


The well-known passivity and small-gain theorems [3] can easily be derived as special cases 
of Theorem 2.2. 

In order to study stability of feedback systems with more general supply rates, 
we state an observability assumption. 

Assumption 1: There exists some T > 0 and a non-negative continuous 

function a : R R, with a(0) = 0 and a(a) > 0 for a > 0, such that for identically zero 
input and any initial state xo , we have 

/ y'(0v(0* ^ a (l! *o||). D 

Jo 

(For finite-dimensional linear time-invariant systems, this is equivalent to the standard def- 
inition of observability.) 
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Theorem 2.3 Suppose that 

(a) A given system G has a supply rate Wi := (y, P\y) + w c (u,y ) + (it, R\u) with 

ft i 

/ w l (u,y)dt>0, 

Jo 

(b) H is a system with supply rate w 2 := (y, P 2 y ) + w c (y,u) + (it, R 2 u ) and 

ft i 

/ w 2 (u,y)dt> 0. 

Jo 

If R 2 + Pi and P 2 + Ri are negative definite , and both G and H satisfy Assumption 1, then 
the origin is an asymptotically stable equilibrium point of the closed loop system (Figure 1 ) 
with zero external inputs. 

Proof: We will demonstrate a Lyapunov function for the closed loop system. Let 
Si and S 2 be storage functions for G and H respectively. Define r = (ri,r 2 ) and similarly y. 
We have iti := r\ — y 2 and u 2 := r 2 + y 2 . It is clear that the interconnected system (Figure 
1) is dissipative with storage function S := S\ + S 2 and supply rate 

w{r,y) = u>i(ui,yi) +w 2 (u 2 ,y 2 ). 

We will show that S is in fact a Lyapunov function for the closed loop system. 

Since 

S(xjo)+ [ w{r,y)dt > S(xi) (5) 

Jo — 

it follows that for ri = r 2 = 0, 

S < iwi(iti,yi) + w 2 (u 2 ,y 2 ) 

= wi(-y 2 ,yi) + w 2 (y u y 2 ) 

= y 2(^-1 + -^ 2 ) 3/2 + y[{P\ + Ri)y\ 

< 0 

since the terms arising from the cross rate w c ' cancel. Thus , S is negative definite (Assump- 
tion 1). 

It now remains only to show that S is a positive definite function. 

S{xo)> S(xi) - [ w(r i y)dt > f -y 2 (Ri + ^ 2 ) 3 / 2 ^ + / -y[(Pi + Ri)yidt. 

Since both plant and controller satisfy Assumption 1, there exists a positive definite function 
/ 3 with 

S(x 0) > fi{\\ ®oll) 

and so S is a Lyapunov function for the closed loop system. The result follows. □ 

Theorems 2.2 and 2.3 will be used to derive controllers for uncertain structures. 
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3 Rate Sensors 


Suppose we have a structural control system (1) with only rate sensors i.e. Cj — C a — 0. 
It is known that if C' v = F then the system transfer function is positive, or equivalently, 
the system is (0, /, 0)-dissipative [1]. The passivity theorem can be used. In particular, any 
strictly positive stable controller will lead to a stable closed loop system. It is not necessary 
to determine the model (1) beyond ensuring that C' v = F and that the dynamics are second 
order. 

In this section we extend this result to more general configurations. 

Theorem 3.1 Consider (1) with Ca = C a = 0 and suppose that there exists an operator 
Q : U i— » Y such that 

C'Q. = F. (6) 

Then, for some p 2 > 0, the system (1) is {—p 2 1 ,Q ,0)— dissipative. 

Proof: For arbitrary p > 0, define the supply rate 
w(u, y ) := ( y , -p 2 y) + (y, Qv-)- 


Define 


t < T 
t > T 


ux(t) := -Pru(t) := | ^ 

and let yr{t) be the output which corresponds to the input ur(0- 
[T fT 

/ w(u, y)dt = / — p 2 {yT,yr) + {yT,Qur)dt 

JO Jo 

*oo too 

= / -p 2 {y T ,yT) + (yT,Qu T )dt + / p 2 (yT,yr)dt 

Jo J T 

poo 

> / -p 2 {yT,yr) + (yT,Qur)dt 

Jo 

1 f 00 

= — / —p 2 (yT,yT) + (yT,Qwr)dv 

Z'K J — oo 

where y denotes the Laplace transform of y. 

Defining v(joj) := — Mio 2 + Djio) -1 Fut(jw), so y = C v v , 


(7) 


f T w G (u,y)dt > J- / -p 2 (C v v,C v v) + {y,Fwr)dw 

Jo Z 7 T J — oo 

= — f°° v(ju)*[-p 2 C v ’C v + —(K - Mu 2 + Djw)]v(jw)dw 
27T J - CO JOJ 

= — T v{]uY[-V 2 CdC d + D)}v(]u)du,. 

2 7 T J- oo 
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If we choose p 2 small enough so that p 2 C v 'C v < D, the expression in square 
brackets is positive definite for all u>. □ 

The following stability theorem is now immediate. 

Theorem 3.2 Assume we have a second order system with rate sensors, and C' V Q = F, as 
in the statement of the previous theorem. If a given controller H is (1) externally L 2 -stable 
and (2) QH is a positive transfer function, then the closed loop system is stable. 

Proof: The system is (— p 2 I, Q, 0)-dissipative for some p > 0 (Theorem 3.1). Let k be the 
controller gain and choose a, b so k < b 2 /a 2 and b 2 < p 2 . We have that the controller is 
(— a 2 /, Q' , 6 2 /)-dissipative. The result now follows from Theorem 2.2. □ 

This result is more general in several ways than previous results depending on 
positivity of a structure. The plant is not constrained to be square and the operator Q is 
not required to bear any relation to a Lyapunov function for the system. However, the rate 
sensors and force actuators do need to be sufficiently “collocated” so that C' V Q = F for some 

Q- 

4 Structures with Displacement Sensors 

In many applications, sensing is performed with displacement sensors. If the control is 
still input as a force, then the system transfer function is by no means positive, even the 
extended sense of (0,Q,0)- dissipativeness discussed in the previous section. We would like 
to obtain a robustness result similar to that obtained in the previous section. One approach 
to this problem is to use multipliers (Figure 2). Figure 2 can be shown to be equivalent to 
Figure 1. The multiplier function L is chosen so that a standard result such as the passivity 
theorem may be used for the transformed systems. Unfortunately, for a given plant, there 
is no guarantee that a suitable multiplier function exists, nor any clue of how to choose 
the multiplier function. Furthermore difficulties associated with ensuring causality of the 
transformed systems may arise. Details can be found in [3]. 

Looking at the situation from a physical viewpoint, replacing rate sensors on a 
structure with displacement sensors does not change the internal dynamics of the structure. 
Intuition tells us that the system should remain dissipative. If suitable supply and storage 
functions do exist, so that the system (1 ) with displacement sensors is dissipative, Theorem 
2.3 can then be applied to design robust controllers. 

In order to keep the discussion simple, we consider the situation where only 
measurements of displacement are made i.e. C v = C a — 0. We note first that the definition 
of supply functions can be extended to include functions which involve derivatives of the 
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inputs and outputs. The input and output spaces have to be defined appropriately so that 
the supply function is well-defined. 

Theorem 4.1 Consider the second-order system (1) with C a — C v = 0, assume that there 
exists Q so C' d Q = F . Let c 2 > 0 be a lower bound on system damping so that c 2 M < D . 
Define the cross rate 

w c (u,y) := c 2 (y,Qu) + (y,Qu). (8) 

Then there is a number p 2 > 0 so that this system is dissipative with respect to the supply 
rate 

w G (u, y) = -p 2 (y, y) + w c (u, y). 


Proof: First note that y exists for all outputs y so that w G is a well-defined supply 
function. As in Theorem 3.1 and in [3] we will show that 

f T 

/ w G (u, y)dt > 0 
Jo 

for arbitrary T > 0 by showing that its transfer function is positive. Dissipativeness of the 
system with this supply rate will then follow from Theorem 2.1. 

Defining ur(t) Pru(t) as in equation (7), let yr(i) be the output which corre- 
sponds to the input ux(t). 
fT fT 

J w G (u,y)dt = I -p (yT,yT) + c 2 (y T ,Qu T ) + (yT,QuT)dt 

fOO roo 

= / — P 2 {yT)VT ) + c 2 (yT, Q u t ) + (yT> Qu?)dt + p (yr, yr)dt 

J 0 v T 

y*oo 

> / —p 2 (yr, 2 /r) + c 2 (y T , Qut) + ( y'r , Qux)dt 

Jo 

f ~P 3 (yf, yr) + c 2 (yr, Qut) + Qur)dw. 

J — oo 


1 t°° 

27T 


Defining v(jlj) := ( K — MuP + Djcj) 

! T w G {u, y )dt > i [°° v{]u)*[-p 2 C d 'C d + c\K - Mu 2 ) + Du 2 )v{ju)du. 

JO 27T J— oo 

If we choose p 2 so that p 2 C d C d < c*K, the expression in square brackets is 
positive definite for all u> and for arbitrary T, 

rT 


f w G (u,y)dt > 0 
Jo 


as required. □ 


8 


Theorem 4.2 Consider the structural system defined in the previous theorem. Assume ob- 
servability. The closed loop system (Figure 1) will be stable if a linear controller is (a) stable 
(b) dissipative with respect to the supply rate 

w(u, y) = w c (y, u) = c 2 (u, Qy ) + (u, Qy) (9) 

where Q,c 2 and w c are as defined in the statement of the previous theorem. The input space 
U is restricted to inputs with u 6 C 2 (U). 

Proof: Define the positive number p 2 as in the previous theorem. Since the controller is 
stable, it is (— a 2 1, 0, 6 2 7)-dissipative where we choose b 2 < p 2 and b 2 /a 2 large than the 
controller gain. Choose an observable realization for the controller. If condition (b) also 
holds, then it follows that the origin is a stable equilibrium point of the closed loop system 
with Ui — u 2 — 0(Theorem 2.3 ). Since both plant and controller are linear and finite 
dimensional, it follows that the closed loop system is also externally stable. □ 

Errors in the model of the system stiffness K or actual damping larger than the 
estimate will not destabilize the closed loop system. Other than ensuring C' d Q = F, the only 
knowledge of the structure required in designing a stabilizing controller is a lower bound on 
the damping so that the controller satisfies (9). If the actual system damping is greater than 
the estimate, the closed loop system will still be stable. 

Note that this result could also have been obtained through defining the multiplier 
function L = Q(c 2 + s). We feel however that there are several advantages to using dissipative 
system theory. The first one is that supply rates such as (u,y) may be defined without 
the problems of causality and stability associated with the multiplier function s. More 
importantly, we think that the the choice of supply rate is clearer when it is known that the 
system “dissipates energy” in some sense. 

5 Example: Second Order Controllers 

In order to illustrate the ideas discussed above, consider the second order dynamic system 
with displacement sensors and C' d Q = F for some Q as discussed in the previous section. 

Suppose we wish to control this system with a combination of direct feedback 
and the displacements of a second order system. Mathematically, the proposed controller 
has the form 

M 2 x 2 {t ) + D 2 x 2 (t ) + K 2 x(t) = F 2 u 2 (t ) (10) 

y 2 (t) = C d2 x 2 {t ) + Eu 2 {t). 
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Since both plant and controller are stable some readers might expect stability of 
the closed loop system to follow if both plant and controller have positive damping. This is 
not true, as the following simple example illustrates. 

Example : x(<) + 8x(i) + 100x(i) = ui(t) x 2 (t) + .001x*2(<) + .lx 2 (i) = u a(0 

u x = ri — x 2 u 2 = r 2 + x 

While for damping d = 8, the closed loop is stable, increasing the plant damping to 12 will 
lead to an unstable closed loop. 

However, by using the results from the previous section we can obtain a controller 
(10) which is unconditionally stable. 

Theorem 5.1 Suppose we have a observable second order plant (1) with C v = C a = 0 
and (2) there exists Q such that C^Q = F . If the controller (10) is implemented so that 
CftQ* = —F 2 and the bias E satisfies QE > F 2 X t F 2 where X > K 2 1 , then the closed loop 

system is stable . 

Proof: From Theorems 4.1 and 4.2 the result will follow if we can show that the 
given system is dissipative with respect to the supply rate (9) for arbitrarily small c > 0. 

Define u(juj) := (— M 2 uP + D 2 ju + -K2) and set E Q = F 2 XF 2 4“ E 2 

where X,E 2 are as yet undetermined. Using the same frequency domain technique as in 
Theorems 3.1 and 4.1, we obtain 

[ T w c (u,y)dt > c 2 f v* ( jw)Z(]u)v(ju) + ut* {ju)G 2 u t (ju ;) du 

Jq J -00 

where 

Z{]u>) = [j K 2 XK 2 - K 2 ] + a > 2 [D 2 /c 2 + M 2 - M 2 XK 2 - K 2 XM 2 + D 2 XD 2 } + w 4 M 2 XM 2 . 

Choose <? small enough that (1) c 2 M < D and (2) the coefficient of u> 2 is positive definite. 
If XK 2 > I, then Z is positive, and the controller is dissipative with supply rate w c for all 
E 2 > 0. Therefore, the closed loop is stable (Theorem 2.3). That is, the closed loop system 
is stable for all E which satisfy 

QE > F' 2 K 2 -'F 2 .□ (11) 

Does this design have any physical meaning? The answer is positive. Consider 
the special case where the controller and the system have the same number of states, and 
furthermore assume that all states are measurable and the actuators are collocated with the 
sensors. 

M£i(i) + Dx^t) + Kx 1 (t) = u a = 7-x (t) - y 2 (t) (12) 
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Vi = ® 1 - 

Adjust the controller gain so F 2 = K 2 : 


M 2 x 2 (t ) + D 2 x 2 (t) + K 2 x 2 (t) = K 2 u 2 ~ K 2 (r 2 (t) + yi(t)) (13) 


V2 = -K 2 x 2 + Eu 2 = AT 2 Z2 + #(t/i + r 2 ) 
Defining x := (x lt x 2 ) we can write the closed loop system as 


M 0 

0 m 2 


®(t) + 


D 0 

0 D 2 


x{t) + 


K + E -K 2 

k 2 k 2 


x(t) = 


I -E 
0 K 2 


L{t)- 


This is a stable second order system as long as the matrix 


K t = 


K + E -K 2 

k 2 k 2 


is positive definite i.e. it represents the stiffness matrix for two “springs” K and K 2 connected 
in series. This is true if E > K 2 . Noting that F 2 — K 2 and Q — I, this inequality is identical 
to the condition (11) arrived at using the dissipativeness criterion. 

If we have scalar systems, the choice E = k 2 leads to a stable closed loop system. 
Setting the reference inputs to zero, 


= -y 2 = k 2 {x 2 - xi) 

The controller in this case reduces to a spring-mass-dashpot system connected in series with 
the system mass as shown in Figure 3. 

For d 2 = 0, the system is just two spring-masses connected in series. As k 2 is 
increased, the stiffness of the system increases. For d 2 > 0, the system is always stable for 
^2 > 0, and increased amounts of energy are dissipated in the dashpot as d 2 is increased. By 
adjusting the control system parameters, ^2,^2,7712 , we can design the closed system to be 
under-, critically or over-damped. 

6 Conclusions 

We have used results in dissipative system theory and Lynpunov’s Second Method to develop 
stability theorems which generalize the passivity and small-gain theorems. Simple conditions 
for robust, model-independent controllers which stabilize non-passive structures were derived. 
Throughout the paper we have used the second-order form of the dynamical equations. Only 
displacement and velocity feedback were studied. However, the extension to accelerometers, 
and to combinations of these three types of sensors is straightforward. 
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Although for displacement feedback, a lower bound on the system damping is 
required, the controller is otherwise independent of the system parameters. It should be 
emphasized that this is a robustness result with respect to structural uncertainty in the 
absence of measurement uncertainty and other contributing factors. 

Control performance is, of course, dependent on the system characteristics. Knowl- 
edge of the system model can always help improve a controller design. Future work will ex- 
plore combining satisfaction of the appropriate dissipative condition with design techniques 
such as LQG regulator and methods. Finally, the controller has been formulated from 
the continuous-time setting. Actual implementation of the controller, however, most likely 
requires usage of a digital computer. In future work, effects of sampling and time delays 
will be addressed. Other practical issue that can also affect the control performance such as 
measurement noises, and actuator and sensor saturation limits will be investigated. 
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Figure 2: Feedback System with Multipliers 
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Figure 3: A simple spring-mass system with a simple spring-mass controller 
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